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Abstract 

We study the Fermi-edge singularity appearing in the current-voltage characteristics for resonant 
tunneling through a localized level at finite temperature. An explicit expression for the current 
at low temperature and near the threshold for the tunneling process is presented which allows to 
coalesce data taken at different temperatures to a single curve. Based on this scaling function for 
the current we analyze experimental data from a GaAs-AlAs-GaAs tunneling device with embedded 
InAs quantum dots obtained at low temperatures in high magnetic fields. 
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The generic response of a fermionic many body system to the sudden appearance of a 
local perturbation is the Fermi-edge singularity (FES). This response can be observed, e.g., 
in experiments probing X-ray absorption in metals^ or resonant tunneling through localized 
levels.-^ In the latter, as depicted in Fig. ^ the Coulomb interaction with a charge in the local 
state |e$) acts as a one-body scattering potential for the electrons in the leads. The change 




FIG. 1: Schematic sketch of the tunneling process through a local state |ej). 

of occupation of the local levels during the tunneling process generates sudden changes 
of this scattering potential leading to characteristic singularities in the I—V curves at the 
corresponding voltage threshold.— ttt&L&2 A theoretical analysis of the low energy response, 
i.e. in the vicinity of these singularities, at zero temperature is possible using bosoniza- 
tion techniques.— In this approach the response functions can be expressed as correlation 
functions of boundary changing operators in an equivalent 1+1 dimensional conformal field 
theory— For the tunneling current in the vicinity of the voltage threshold Vo this gives a 
power law, 

I(V, T = 0) oc (D/e(V - Vo)) 1 6{V - V ) , (1) 

with a characteristic high energy cutoff D of the order of the band width. The critical 
exponent 7 at this threshold is completely determined by the scattering potential and can 
be expressed in terms of the scattering phase shifts in the channels coupled to the scatterer 
taken at the associated Fermi momenta . 2l3ll ° The appearance of power laws for the response 
of the system in this approach is a consequence of the absence of an intrinsic scale in the 
low-energy theory for this problem. 

There are several effects, however, which affect this singularity and therefore obstruct 
their direct observation: already at zero temperature the finite lifetime of an electron in 
the local state leads to an intrinsic broadening of this level. Since its energy e,, measured 
relative to the Fermi energy, depends on the voltage as ej oc e(V — Vo) this broadening 
directly modifies Eq. (JTJ). 3 In addition, nonequilibrium effects causing dissipation have to 
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be taken into account for a complete description of the system with a voltage bias across 
the tunneling barrier. This implies further smoothing of the FES due to the resulting finite 
lifetime of states in the leads , ^i 1 ?! 14 Even without dissipation the response functions such as 
/(V) of the system are exponentially damped at finite temperature, i.e. away from criticality. 
However, as the temperature is lowered, Eq. (JJJ has to be approached, which is reflected by 
a divergence of the maximum of the measured current as I(V max , T) oc T -7 , again limited by 
the intrinsic broadening of the local level. Here, the critical exponent 7 is the same as the one 
governing bias dependence of the T = response. Finally, the derivation of the expression 
for the tunneling current in the bosonization approach relies on a linear dispersion of 
the lead electrons. While this approximation is valid close to the threshold, nonlinearities 
or fluctuations of the density of states (DOS) due to the presence of impurities in the leads 
will modify the current for finite detuning AV = V — Vq of the voltage. These corrections 
are a property of the probe and independent of the temperature. Therefore, they can be 
distinguished from the FES if its dependence on both temperature and voltage, at least at 
low energies, is known. 

Neglecting other than thermal broadening a closed expression for this dependence can 
be derived for lead electrons with linear dispersion: at finite temperature T and at a finite 
detuning AV of the voltage from the threshold the response of the system will depend on the 
ratio eAV/ksT of these energy scales reproducing the power laws mentioned above for zero 
temperature and at the maximum, respectively. For non-interacting electrons the extension 
of Eq. (^) for the X-ray absorption problem to finite temperatures can be achieved by taking 
into account the Fermi distribution for the initial state o 1 ^ 1 ^ (see also Ref.ll7l). Alternatively 
it can be obtained by incorporating temperature in the bosonization procedure of Ref. 
Within the latter approach the analysis is easily generalized to the case of interacting ID 
fermions described by the Luttinger Hamiltonia n 11 ! 18 : the Luttinger parameter K p only 
enters the exponent 7 in the final expressions Q and To be specific, we consider a 
single channel of non-interacting fermions in the emitter of the tunneling device coupled to 
a local level at the origin. The problem can be formulated as a one-dimensional scattering 
problem in terms of a left moving (chiral) fermion field ip o n the interval [-L, L] with Fermi 
velocity vp described by the Hamiltonian (see e.g. Ref. I 111) 

H= f dx Uv^-^j) + S(x)Uip^blbi\ + e-blbi (2) 



/ely, 
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(bi annihilates the electron in the local state). Using Fermi's golden rule with the tunneling 
operator (^(O)^ + h.c.) the current can be expressed in terms of the Green's function 

Gf(t) = ((WV>(0,t)) + WV(0,0)) T • (3) 

Here, 1A{t) is the unitary operator describing the change of the boundary condition at x = 
due to the tunneling of an electron. It relates the Hamiltonian (j2j) in the sectors with the 
local level occupied or empty by means of a canonical transformation.^ Representing the 
fermion ip oc exp(— i\f^n<\)) in terms of a left-moving boson (f>(x, t) with mode expansion (the 
zero mode O does not contribute to (J3J)) 

oo 

<Kx,t) = 0o + V -77= {a m e- imn{x+VFt)/L + h.c.) 

m=l v 

the Green's function is (5 = U/2vf is the Fermi phase shift) 

{/ r \ 2 oo 

( 1 + - ) J2 - I 2 Ra m ) T (cos (rrnrt/L) - 1) + e " im7r '/ L - l] 

For T = the ground state expectation value of the bosonic occupation numbers vanishes, 
( a m a m)o — 0, and this expression leads directly to ([T)) . 10 ^ 1 At finite temperatures the Bose- 
Einstein distribution of the bosonic occupation numbers (a^a^x is used instead.J^ Within 
this approach we obtain (7 = —25/n — (5/n) 2 ) 



00 
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I(AV,T) oc / dte ie£XVt G(t) 



- /, - , u m (4) 



B — s ^ '7 



— Im 

for the low temperature response (see also Refs. Here (3 = 1/ksT and B is the beta 

function. 

This expression reproduces the power laws mentioned above: (i) at voltages sufficiently 
above the threshold, for eAV ^> fc^T, the power law singularity Eq. (P) emerges from 
Eq. (ii) For positive 7, the tunneling current becomes maximal at some finite value of 
(3eAV which approaches for 7 — > 1 and diverges as 7 — > 0. Since a given FES is described 
by a T-independent exponent 7, this implies that the position of the maximum in the I-V 
curve moves to higher voltage with the temperature Vm ax — Vq oc T. At this maximum, the 
only temperature dependence of the current is in the pre-factor, i.e. J max oc T~ 7 . (iii) Finally, 
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for 7 — > 0, i.e. no scattering of the electrons off the local level irrespective of its occupation, 
Eq. (jlj) reduces to a thermally broadened current step I(AV,T) oc 1/ (1 + exp (— /3eAV)). 

In the following we apply this result to the analysis of edge singularities observed in 
tunneling experiments through localized levels in InAs quantum dots in a strong magnetic 
field. These self-assembled InAs quantum dots are sandwiched between two AlAs tunneling 
barriers. They have a height of approximately 3nm and a diameter of 2r ~ 10 — 15 nm£. 
The emitter and collector consist of a 15 nm thick undoped GaAs spacer layer followed by a 
n-doped GaAs buffer with graded doping: a 10 nm thick layer of n~ doping (1 • 10 16 cm -3 ), a 
10 nm thick layer of n doping (l-10 17 cm~ 3 ) and a 1 jum thick layer of n + doping (2-10 18 cm -3 ). 
We have studied two samples with nearly the same growth parameters; only the tunneling 
barrier thicknesses are different. One has AlAs tunneling barrier thicknesses of 4 and 3nm 
and we show here measurements at a magnetic field of 14.9 T. The other sample with more 
asymmetric barriers (5 and 2 nm) was studied at magnetic fields up to 28 T. In this paper we 
look at the tunneling direction where the electrons tunnel first through the thicker barrier 
onto the dot and leave it through the thinner barrier. 

The thickness of the spacer layer (15 nm) of the investigated samples was carefully chosen 
to allow for a description of the emitter as a three dimensional electron gas even at the 
presented bias voltages, different from the structures used in Ref. £l, A comparison of 
the I-V characteristics to a test sample with a much thicker spacer layer (100 nm) showed 
distinct differences. For this test sample strong peaks in the I-V characteristic due to the 
formation of 2D subbands in the emitter right before the tunneling barriers evolved. In 
Ref. such a 2D emitter behavior was also observed for the same 100 nm thickness of the 
spacer layer. The absence of such peaks in all our samples with a 15 nm spacer layer confirms 
our assumption of a 3D emitter. Additional support fro this assumption has been given by 
a detailed analysis of a sample of the wafer with the 4 and 3 nm barriers in Ref. |22j. Both 
the current and the measured shot noise showed excellent agreement with the theoretical 
models for 3D-0D-3D tunneling. 

Expression for the tunneling current implies that the I-V data taken at different 
temperatures T - after rescaling the current as IT 1 and the voltage as eAV/ksT - should 
lie on a single scaling curve. Here, the only free parameter should be the edge exponent 
7. In practice, however, there are various additional effects to be considered: first, only a 
fraction a of the total applied voltage drop occurs between the emitter and the quantum 
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dot, i.e. 6j = aeAV. The energy-to-voltage conversion factor a can be determined from the 
thermal broadening of the FES. Second, as already mentioned in the introduction, the local 
level in the quantum dot has an intrinsic width Tj due to its hybridization with the states in 
the emitter. The effect of these mechanisms on the broadening of the I-V curve is different 
and non-equilibrium condition o 1 ?' 14 have to be taken into account for the latter. Here we 
include both by introducing an effective temperature k B T' = a/ (k B T) 2 + Tf in Eq. (JU). 
Based on the theoretical expression we have analyzed data taken with the first sample in a 
magnetic field of 14.9 T at various temperatures between 0.3 and 1.0 K (Fig. 12(a)). Due to 
the Zeeman-splitting of the local level in the dot the FES appear in pairs (see also Refs.lfif^L 
Hence the data have to be described by the sum of two contributions (jlj) with different edge 
exponents 7j, 7^ but the same conversion factor a and intrinsic broadening IV Unlike in 
bulk InAs the Lande factor of the quantum dots is positive.—^ As a consequence we identify 
the FES at lower bias voltage with tunneling of spin-j electrons. Performing a fit to the 
experimental data we obtain 7j = 0.27, 7^ = 0.46, a = 0.15 and Tj = 0.36 K. Now, the 
FES from tunneling of electrons with spin j (corresponding to the peak in the I-V curves 
at lower bias) can be isolated by subtracting the theoretical contribution from the other 
spin projection. Rescaling the axes as described above the I-V data are indeed seen to 
collapse onto a single curve given by (@J) (see Fig. 12(b)). The scatter observed at higher 
voltages is caused by fluctuations in the local DOS in the emitter. As argued above these 
fluctuations should lead to temperature independent deviations of the experimental data 
from the theoretical prediction. This behaviour is shown in Fig. El where the raw I-V data 
(i.e. without rescaling of the axes) near the Zeeman-split FES at two different temperatures 
are shown. Taking the ratio of the experimental data to the theoretical fit the temperature 
independence of the fluctuations in the current becomes manifest, see Fig. E(b). Although 
the data show some scatter which is most likely caused by noise in the experimental data, it 
is astonishing that in all five experimental curves the same small deviation in the theoretical 
curve can be seen. 

As has been reported before,^ the edge exponents characterizing the FES in the tunneling 
current depend strongly on the magnetic field. To understand this in terms of the parameters 
describing the dot and its interaction with the emitter we have to go beyond the generic 
low-energy Hamiltonian (J2J): the broadening of the localized level and the Fermi phase 
shifts determining the edge exponents in (JTJ), (jlj) need to be computed within a suitable 



6 



235 



240 245 

V [mV] 



250 




FIG. 2: (a) I-V characteristics exhibiting the spin-split Fermi edge singularity for a magnetic field 
of 14.9 T at temperatures T = 0.3 K, 0.35 K, 0.45 if, 0.55 K, 0.6 K, 0.75 K, 1.0 K (for clarity a 
T-dependent offset is added to I), (b) Same data for the FES at V ~ 238ml/ with current rescaled 
as IT' 1 vs. voltage in units of the effective temperature T 1 (see main text) for a collapse onto a 
single curve. The full line is the best fit of Eq. (JIJ to the data. 

microscopic model. As discussed above the emitter can be described as a 3D electron gas 



in the half space z < (see, e.g., Ref |24j for the description of an effective 2D electron 
system). In the experiment the transverse motion of the electrons is quantized by the 
magnetic field B\\z. Only the lowest Landau band (LLB), n = 0, is occupied and the 
electrons populate one-dimensional channels of given angular momentum m > with Fermi 
momentum ftkp a for spin projection o =|, j.— The radial wave functions in these channels, 
ipn=o,m(p, 0) oc p m exp(— zm0 — p 2 /4£g), depend on the magnetic field through the magnetic 



length £q = \Jh/eB. In the experiments, £o is comparable to the lateral size r$ of the 
quantum dot. 

Assuming a Gaussian wave function for the electron on the isolated quantum dot the 
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FIG. 3: (a) I — V data for a Zeeman-split edge singulartiy at B = 14.9 T, T = 0.3 K (upper curve) 
and T = 1.0 K (lower curve). The dots show the experimental data, the solid line is the theoretical 
fit. (b) Ratio of experimental data and theoretical fit for temperatures T = 0.3 K, 0.45 K, 0.6 K, 
0.75 K, 1.0 -PC showing the (essentially T-independent) fluctuations in the local DOS of the emitter. 

intrinsic broadening Ti can be computed in perturbation theory from its overlap with the 
single electron states in the LLB of the leads (in the present geometry the broadening will 
be dominated by the overlap with the collector states). The field dependence of the wave 
functions ipo m leads to a linear growth of the broadening with the magnetic field with a 
change in slope around £q ~ r . This agrees with the qualitative behaviour of Tj as obtained 
from our fits to the experimental I-V data: analyzing the evolution of the FES with the 
field we find that Tj varies between 0.36 K at B = 14.9 T and 1.5 K at B = 28 T. as a 
consequence, the corresponding lifetimes of particles in the local state are well above the 
thresholds for the observation of the FES.^. 

Finally, we want to discuss the magnetic field dependence of the edge exponents 7 CT 
characterizing the Zeeman-split FES. In the tunneling experiment the scattering potential 
of the dot affects several channels labelled by the angular momentum m and spin a of the 
electrons. Alternation of the occupation of the quantum dot implies changes S m (k) in the 
Fermi phase shifts experienced by the electrons in channel m. As in the simplified model 
(J2J the edge exponents can be expressed in terms of these changes (angular momentum 
conservation implies that only m = electrons contribute to the tunneling matrix element)^ 
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7<x = 




(5) 



The scattering potential for the electrons in the emitter is composed from the band bending 
due to the presence of the semiconductor-insulator interface at z = and the screened 
Coulomb potential of a charge on the quantum dot - if it is occupied. In the ID Landau 
channels with given angular momentum m both can be described by effective potentials 
ve KZ decaying exponentially into the emitter with range k _1 being the Debye radius. The 
resulting phase shift picked up by the ID electrons with momentum Hk perpendicular to the 
boundary is 



where o-Fi is a hypergeometrical function. For the empty dot the effective potential v is the 
band bending Vb depending on the doping profile. For the charged dot v — v & + iVm where 
v n>m is obtained by projecting the screened Coulomb potential V(p, z) (see e.g. Ref. |3J) to the 
radial wave functions if) n ,m{p, 4>) of the electrons in the n th Landau level^. Both the effective 
potential obtained by this projection and the Fermi momenta kp a in the Landau quantized 
channels of the emitter are functions of the applied magnetic field. Through Eq. (JHJ) this 
leads to the strong B dependence of the edge exponents -y a observed in experiments with 
3D emitter. 

In Ref. Q the FES was analyzed based on the linear behaviour of © for small k. In this 
approximation the difference of the exponents 7j — 7j_ is proportional to the corresponding 
Fermi momenta which can be determined from the ID DOS in the Landau bands (broadened 
by disorder). As a consequence, the difference of the edge exponents should grow mono- 
tonically with B. This behaviour is indeed observed in the experiments at sufficiently high 
magnetic fields B > B c w 20 T. Below B c , however, the experimental data presented in 
Ref. Q indicate that 7-f — 7j_ changes sign, see also Fig. HJ According to Eq. (J3J) such a cross- 
ing requires that 5(kpi) = S(kpi) for B c . Clearly, such a feature requires a non-monotonic 
^-dependence of the Fermi phase shift and can only be captured when the full nonlinear 
expression (JHJ) for the scattering phase is taken into account. In Fig.HJ we present a theoret- 
ical prediction based on this expression. Both Landau level broadening (affecting the field 
dependence of the Fermi momenta) and bending of the Landau band in the emitter at z = 
have to be taken into account. The latter is found to be close to the threshold of formation 





9 




FIG. 4: Field dependence of the exponent 7 for both spin polarisations. For B > 20 T fitting of 
Eq. to the FES for the minority spin is hampered by its decreasing weight. This results in the 
large errors displayed for 7j_. Full lines are theoretical predictions based on the phase shift © of 
the exponentially screened scattering potential. 

of a single electron bound state at the interface. This is consistent with our picture that at 
very high bias voltages a potential well is formed in the emitter close to the barriers. At 
the same time this limits the validity of the simple description for the scattering potential 
in the emitter in terms of a decaying exponential used in the derivation of Eq. ©• Within 
this simple model the qualitative S-dependence of the experimental data near B c is repro- 
duced. The splitting of the edge exponents at B < B c , however, is underestimated by the 
prediction. Note that the z-dependence of the scattering potential for the electrons in the 
emitter can, in principle, be determined from the field dependence of the edge exponents. 
For this, however, they have to be extracted from the experimental data with sufficient 
accuracy. Here the quality of data is limited in particular by the low weight of the FES for 
the minority spin j. at fields B > B c . 

In this Paper we have analyzed FES observed in resonant tunneling experiments through 
InAs quantum dots between 3D leads in strong magnetic fields parallel to the current using 
both its temperature and bias-voltage dependence. Based on an explicit expression (JU) for 
the I-V characteristic of the FES data taken at different temperatures have been collapsed 
into a single scaling form. From these data the effect of the temperature-independent fluc- 
tuations in the local density of states in the emitter on the current has been identified. 
Furthermore, the parameters describing the interaction between the dot and the lead could 
be extracted. Based on a microscopic model for the tunneling device the magnetic field 
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dependence of the broadening of the local level on the quantum dot has been described. In 
addition the exponent characterizing the FES for both spin projections has been computed 
taking into account the non-linear momentum dependence of the Fermi phase shift (JBJ), re- 
producing the qualitative behaviour found in the experiment, in particular the change in 
sign of 7f — 7j_ at a magnetic field B w 20 T. For an improved description, both the effect of 
non-equilibrium condition s 13 ^ 4 has to be considered and fine-tuning the scattering potential 
for the electrons in the emitter used to compute the Fermi phase shift is necessary. The 
latter would amount to introducing introduction of additional fitting parameters, though. 
On the other hand, we note that, this potential can be reconstructed in principle from 
the experimental data provided that the edge exponents can be extracted with sufficient 
accuracy. 
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